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ABSTRACT 

This  report  concerns  two  so-called  inverse  problems  of  mathematical 
physics.  These  are:  (i)  the  problem  of  determining  a second-order  differ- 
ential operator  (in  a normal  form)  on  the  half-axis  from  its  spectral  function; 
and,  (ii)  the  problem  of  determining  a hyperbolic  boundary  value  problem  of  a 
special  form  in  a (non-characteristic)  half-plane  from  its  response  on  the 
boundary  to  a unit  impulse  at  some  reference  time  t = 0 (boundary  value  of 
the  Riemann  function). 

We  solve  problem  (ii)  by  a natural  approach,  and  then  indicate  how  the 
solution  of  problem  (i)  follows  from  the  solution  of  problem  (ii).  Our  solu- 
tion of  problem  (ii)  is  constructive,  and  we  obtain  stability  of  the  solution 
under  perturbation  of  the  data,  in  a well-defined  sense. 

For  problem  (i),  we  obtain  the  well-known  result  of  Gel 'f and  and  Levitan, 
in  the  sharp  formulation  given  by  Levitan  and  Gasymov  ( [6] ) . 
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Significance  and  Explanation 


This  report  concerns  an  inverse  problem  of  wave  propagation.  In  general, 
inverse  problems  for  dynamical  systems  require  that  the  dynamical  laws  of  a 
system  be  recovered  from  information  about  the  behaviour  of  the  system,  such 
as  its  response  to  a specific  stimulus,  or  its  long-time  asymptotic  behaviour. 
Direct  problems  of  dynamical  systems,  which  require  that  the  behaviour  of  a 
system  be  calculated  from  the  dynamical  laws , assumed  known , have  a long 
history,  and  comprise  the  bulk  of  classical  mathematical  physics.  Of  course 
our  understanding  of  physical  systems  proceeds  by  the  solution  of  inverse 
problems:  we  observe  the  behaviour  of  a system,  and  try  to  infer  its  dynamics. 
Nevertheless,  such  inverse  problems  are  usually  approached  by  a sequence  of 
guesses  about  the  dynamical  laws  perhaps  with  some  free  parameters  followed 
by  solutions  of  the  corresponding  direct  problems,  and  attempts  to  adjust 
the  free  parameters  so  that  the  predicted  behaviour  of  the  system  matches  its 
observed  behaviour.  Frontal  attacks  on  specific  inverse  problems,  on  the 
other  hand,  have  been  relatively  rare. 

In  this  report,  we  consider  a very  simple  model  inverse  problem  of  wave 
propagation  in  a one-dimensional  inhomogeneous  medium.  Examples  of  physical 
problems  to  which  our  methods  apply  are:  (i)  the  description  of  the 
(inhomogeneous)  density  distribution  of  a vibrating  string  from  observation 
of  the  motion  of  one  point  along  its  length;  (ii)  synthesis  of  a cable  of  a 
restricted  type,  having  prescribed  transmission  characteristics  for  signals 
of  various  frequencies. 

In  actual  applications , stability  is  an  important  attribute  of  any  compu- 
tation scheme.  One  wants  to  ensure  that  small  errors  in  the  input  data  do  not 
result  in  huge,  uncontrollable  errors  in  the  solution;  otherwise,  the  scheme  is 
generally  useless.  We  show  that  our  reconstruction  of  the  dynamical  law  of 
our  simple  model  system  is  stable,  in  a precise  sense,  against  small  perturba- 
tions of  the  "observed"  data;  in  fact,  we  give  an  approximation  scheme  with 
explicit  error  bounds. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


INVERSE  BOUNDARY  VALUE  PROBLEMS 


* AND  A THEOREM  OF  GEL'FAND  AND  LEVITAN 

W.  Symes 

SI.  Introduction. 

This  report  concerns  two  so-called  inverse  problems  of  mathematical  physics.  These  are: 
(i)  the  problem  of  determining  a second-order  differential  operator  (in  a normal  form)  on 
the  half-axis  from  its  spectral  function;  and,  (ii)  the  problem  of  determining  a hyperbolic 
boundary  value  problem  of  a special  form  in  a (non-characteristic ) half-plane  from  its  re- 
sponse on  the  boundary  to  a unit  impulse  at  some  reference  time  t = 0 (boundary  value  of  the 

' 

Riemann  function). 

We  solve  problem  (ii)  by  a natural  approach,  and  then  indicate  how  the  solution  of 
problem  (i)  follows  from  the  solution  of  problem  (ii).  Our  solution  of  problem  (ii)  is  con- 
structive, and  we  obtain  stability  of  the  solution  under  perturbation  of  the  data,  in  a well- 
defined  sense. 

For  problem  (i) , we  obtain  the  well-known  result  of  Gel'fand  and  Levitan,  in  the  sharp 
formulation  given  by  Levitan  and  Gasymov  ( (6) ) . 


* 


§2.  Notation  and  Statement  of  Results. 

+ _ k ** 

We  will  write  P = [0,«>)  c ]R  throughout.  We  use  the  common  notations  C (U),  C (U)  , 

CQ  (U)  , 8^  (U),  S'  (u)  , D'  (U)  for  the  space  of  k-times,  respectively  infinitely  differentiable 

functions,  respectively  those  with  compact  support,  with  the  usual  Freche’t  topologies,  and 

their  strong  duals.  Here  U c p n may  be  closed  or  open. 

We  denote  by  (U)  the  collection  of  functions  in  Cm  1 (U)  whose  m—  partial 

derivatives,  which  a priori  exist  as  distributions,  may  be  identified  with  locally  absolutely 

integrable  functions  on  U c pn  , i.e.  lie  in  L?'  (U) . We  give  Wm,^(U)  its  usual  Freche't 

loc  loc 

topology . 

If  U c p , denote  by  y}’'  ^(U*V)  the  space  of  continuous  functions  on  U x V whose 
loc 

(distributional)  partial  derivative  in  the  first  (second)  variable  may  be  identified  with  a 

continuous  function  of  the  second  (first)  variable  with  values  in  L*  (U)  (L^  (V)). 

loc  loc 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024  and  the  National 
Science  Foundation  under  Grant  No.  MCS75-17385  A01. 


Recording  to  Fubini's  theorem, 


fc'J'^tOxV)  c W^^tOxV)  . 
loc  Xoc 


11  2 
It  is  clear  how  to  define  U , * (Q)  for  an  arbitrary  open  set  Q c » , since  the 

loc 

definition  is  local,  and  any  p t Q has  a product  neighborhood.  On  the  other  hand, 

suppose  Q is  closed  with  smooth  boundary,  and  let  f e C°  (Q)  n fc'J^Unt  Q) . Let 

2 

p t 3 Q,  and  select  U,V  c jr  open  so  that  uxv  is  a neighborhood  of  p in  H . Ihen 
for  each  x e U (y  e V)  , extend  the  partial  derivative  D2f(x,0  (D^ff-.y))  to  a distri- 
bution on  V (U)  by  requiring 

< D2f  (x,  •),♦>  = ( Djf  (x,  •)  ;$|v  > for  4>  e c”(V) 

x 

« Dj^ft-.y),*)  = < Djf  (*,y)  ,$1^  > for  4>  e c”(U) 

y 

where  V = (ye  V:  (x,y)  e Q}  (U  = {x  e U : (x,y)  e g}). 

X y 

Vte  declare  that  f e W ^^(Q)  if  and  only  if,  for  each  choice  of  p,U,V  as  above,  each 

x c U (ye  V),  the  partial  derivative  D2f(x, •)  (D^ff'.y))  may  be  identified  with  a 

locally  integrable  function  on  V (U) , and  the  map  x t*  ^ffx,*)  (y  D^ft-.y))  is 

continuous,  i.e.  lies  in  C°  (U;  L^  (V))  (C°(V;  L?'  (0)). 

loc  loc 

The  topology  on  U (Q)  is  given  by  the  C°  (Q) -seminorm  and  the  local  norms  of 


the  form 


sup  / dy  | D f (x,y) | , sup  / dx|Df(x,y)| 

xeK  L y«L  K 


with  K,  L c r compact. 

Finally,  ft  (Q)  < ® > 1(  if  and  only  if  f e C™*”1  (Q)  and  all  m-1—  order 

partial  derivatives  of  f are  in  The  topology  on  W ™ 1 (Q)  is  defined  in  a 

loc  loc 

similar  way  to  that  of  V^'^(g).  We  note  that  WI'’'^(Q)  c W™'^(Q),  and  that  the  topology 

XOC  xoc  xoc 

on  W®'*  is  stronger  than  that  of  W™'1  . 
loc  loc 

Functions  in  (g)  may  be  constructed  in  the  following  way:  suppose  f e ( JO . 

Then,  as  is  easily  verified,  the  function  F : Q •*  I defined  by 

F(x,y)  = f (ax  + by) 

for  some  (a,b)  e *2,  lies  in  w”'1(Q),  for  any  closed  g with  piecewise  smooth  boundary. 

loc 

In  fact,  the  spaces  V®^  enter  the  theory  developed  here  in  precisely  this  way. 

We  shall  use  various  common  notations  for  derivatives  and  partial  derivatives,  such 


as  primes,  subscripts,  D ,,  etc.,  without  comment. 

V Jx 

Our  results  provide  an  alternate  route  for  part  of  the  proof  of  the  following  theorem, 
which  is  a sharp  version,  due  to  Levitan  and  Gasymov  (16]),  of  the  celebrated  theorem  of 
Gel'fand  and  Levitan  ([4]). 

Theorem  I.  A nondecreasing  function  p:F  ->  R is  the  spectral  function  of  a 
boundary  value  problem  on  R + = [0,"°): 

-y"  + (q  - X )y  = 0 
y'(0)+hey(0)=0,  h e R 

with  q W™  (R+),  if  and  only  if  p satisfies  the  conditions: 
loc 

(i)  the  integrals 

/N  do  (A)  cos  f.  x = I^fx) 

converge  absolutely  to  functions  in  ( F+  ) , whe  re 


o(A)  = 


■{ 


p (A) Y\ 

71 

p (X) 


( 2.  o 

A < 0 


m+1  + % 

Moreover,  (lN  : N c S)  converges  as  N -*•  00  in  w^oc  (R  ) to  a function  f 


with 


: (0) 


h . 


(ii)  Suppose  u « L2  (F+  ) has  compact  support.  Let 

'V,  fOD  ~ 

u(A)  = J dx  u(x)  cos  /A  x . 

2 ° 

Then  d t L (F  ;dp),  and 

ro  , 'v*  1 2 

dp  ( A)  | u ( A)  | =0 

• 00 

if  and  only  if  u ; 0 a.e. 

Remark.  In  (ii)  above  and  for  the  rest  of  this  paper,  "cos  /A  x"  denotes  the  entire 

function  of  A whose  value  for  A > 0 is  cos  /a  x . In  particular,  for  A < 0 , 
cos  /a  x = cosh  /|  A |x. 

The  bulk  of  this  paper  is  devoted  to  proving  the  following  two  theorems,  from  which 
Theorem  I follows,  and  which  are  of  interest  in  their  own  right: 

Theorem  II.  An  even  function  h : R -♦  R is  the  boundary  value  of  the  Riemann 
function  : 

hit)  = R (0 , t ,0 ,0 ) , t * 0 
for  a hyperbolic  boundary  value  problem: 


-3- 


u - u + q u = 0 

tt  XX 

u <0,t)  + h u (0  ,t)  HO,  t e B 

X 

with  q e wf  (B+)  if  and  only  if 
loc 

(i)  It  e wT1  <*>  * 2(0)  “ h 

loc 

(ii)  the  kernel  f(s,t)  = j(l*(s+t)  + ^(s-t)>  satisfies  the  condition: 

2 

for  any  T > 0 , there  exists  e(T)  >0  so  that  for  all  u £ L (10,T)), 

/T  ds  /T  dt  u(s)  u(t)f(s,t)  + /T  |u  | 2 >_  £ (T)  /T  | u | 2 . 

0 0 0 0 

Remark . An  easy  compactness  argument  shows  that  the  condition  (ii)  above  is  equivalent 

to  the  assertion  that 

f°  ds  / dt  u(s)  u(t)  f(s,t)  + f°  | u | 2 > 0 

0 2 °+  0 
for  all  u e L (B  ) with  compact  support. 

Theorem  III  (i)  The  collection  of  even  functions  f*  £ wT+^'*(B)  defined 
loc 

by  condition  (ii)  in  the  statement  of  Theorem  II  forms  an  open  set  in  wf+1'1,(B) 

loc  even 

(ii)  for  any  even  f in  wJ^'^dO  ' ^*lere  are  precisely  one  h £ B and  one 

q £ W™'^(B+)  so  that  f(t)  = R(0,t,0,0),  t ^ 0,  for  the  Riemann  function  R of 
loc 

the  boundary  value  problem 

u.  . - u + q u = 0 

tt  XX 

ux<0,t)  + h u(0,t)  = 0 

(iii)  The  map  £ )»  (h,g)  whose  existence  is  implicit  in  statement  (ii),  is  a 

continuous  map  from  the  open  set  described  in  statement  (i)  to  B»W?’*(B+)  . 

loc 

Remark . Theorem  III  is  a uniqueness  and  stability  theorem.  Theorem  I is , of  course, 
more-or-less  well  known,  and  Theorem  II  could  be  deduced  from  Theorem  I.  Our  method  of 
proof,  hcwever,  proceeds  by  means  of  an  iteration  scheme,  with  error  bounds  given  explicitly 
in  terms  of  the  nunbers  e(T)  mentioned  in  Theorem  II,  and  various  norms  of  ? . In 


particular,  we  obtain  the  stability  statement  of  Theorem  III,  which  seems  to  be  new. 


(3.1) 


§3.  Heuristic  Discussion  of  Results 

Problem  (i)  is  concerned  with  the  spectral  function  of  a boundary  value  problem 
-y"  (x)  + (q(x)  - A)y(x)  = 0,  X ^ 0,  \en 
y1  (0)  + h y (0)  =0 

where  h is  some  real  number.  Let  $(x,A)  be  the  solution  to  (1)  selected  by  the  initial 
condition 

<M0,A)  - 1 . 

Ihe  spectral  function  of  (1)  is  a nondecreasing  function  whose  associated  Stieltjes  measure 
properly  weights  the  "eigenfunctions"  $(*,A)  in  the  spectral  resolution  of  the  identity 
for  (1),  which  is  concisely  written 

S(x-y)  = dp  (A)  <Mx,A)  <J>  Cy , A ) . (3.2) 

• 00 

In  this  heuristic  discussion,  we  shall  not  worry  about  making  precise  sense  of  diver- 
gent integrals  such  as  (2);  that  is  done,  in  any  case,  in  standard  textbooks  on  spectral 
theory  of  ordinary  differential  operators,  e.g.  [2],  Ch.  XIII.  Nor  shall  we  make  precise 
smoothness  assumptions  on  q . 

Problem  (i)  is:  given  p , find  the  differential  equation  (that  is,  q)  and  the 
boundary  condition  (that  is,  h),  which  give  rise  to  p . Of  course,  this  involves  de- 
scribing those  nondecreasing  p that  arise  as  spectral  functions  of  problems  of  type  (1). 

Problem  (i)  is  a refined  version  of  the  inverse  eigenvalue  problem  : to  construct 
a differential  operator  of  some  special  type,  cum  boundary  conditions,  having  a given 
spectrum.  This  problem  admits  a leurge  amount  of  non-uniqueness  in  its  solution.  Since 
the  points  of  increase  of  the  spectral  function  of  (1)  exactly  amount  to  its  spectrum, 
a solution  to  problem  (i)  certainly  solves  the  inverse  eigenvalue  problem.  The  spectral 
function  also  carries  normalization  information,  however,  and  this  additional  information 
makes  its  solution  unique. 

For  a history  of  these  problems  consult  [3J  and  references  cited  therein. 

Problem  (ii)  is  concerned  with  a hyperbolic  boundary  value  problem 


92  32 

— 2 j-  + q (x) ) 

3t  3x 


3u 


u (x,t) 


(x,t)  c » X R 


(3.3) 


jjj-  <0,t)  ♦ h u (0 ,t ) = 0 , t e * 


-5- 


The  Riemann  Function  R(x,t;x0,tQ)  is  the  solution  of  the  mixed  problem  obtained  by 
adding  to  (3.3)  the  initial  conditions 

u(x,tQ)  « 6 (x-xQ) 

It  <-V  ■ ° • 

Then  u(x,t)  = R(x,t :x0,tQ) . Again,  q is  "smooth  enough”,  and  we  do  not  worry  for  the 
moment  about  the  sense  in  which  a distribution  satisfies  a mixed  problem  of  this  sort. 

In  fact,  as  explained  in  [5],  R is  a special  distribution  with  well-defined  restric- 
tions to  vertical  (and  horizontal)  lines  such  as  {x=0}. 

We  also  note  that  our  Riemann  Function  is  a derivative  of  the  object  usually  called 
by  that  name;  see  (1),  Ch.  VI,  §15. 

Note  that  the  solution  of  the  inhomogeneous  mixed  problem 

,2 


Il32 


3t 


3u 


3x 


+ q(x) 


u (x,t ) = g (x,t ) , t > t , x 


3x 


(0,t)  + h u(0,t)  5 0 


(3.4) 


u(x,tQ)  = |j-  (u,tQ)  5 0 , x > 0 


can  be  represented  in  the  form 


u(x,t)  = y*  da  /°  dT  f°  dy  R(x,t ,y , r)g (y , v) 


(3.5) 


0 0 

(Duhamel's  integral).  The  system  (3.4)  models  certain  processes  (e.g.  nonuniform  trans- 
mission lines)  subject  to  an  imposed  force  g , in  which  signals,  represented  by  u , 
propagate  with  unit  speed,  and  some  boundary  conditions  are  imposed  at  the  "surface" 
x * 0 . Then  (3.5)  represents  the  "forward"  (t  > tQ)  response  of  the  system  to  the 
impulse  g . If  we  choose  in  particular  g(x,t)  - 6 (x-xQ ,t-tQ) , then 


u(x,t)  - J*  R(x,t,x  ,o  )dc,  t > t 


0'  0' 

t0 

we  conclude  that 

ds  £(s)  5 f ds  R(0, s,0,0), 

0 0 

represents  the  response  at  time  t > 0 at  the  surface  x - 0,  to  a unit  impulse  applied 


t > 0 


at  t « 0,  x 


0. 


-6- 


Problem  (ii)  is:  given  f , find  the  boundary  value  problem  (3.3),  that  is,  find 
the  function  q and  the  number  h . Otherwise  put,  we  are  to  recover  the  dynamics  of 
the  system  from  a knowledge  of  its  response  along  the  "surface"  x = 0 to  a unit  im- 
pulse, also  applied  at  the  surface. 

Problem  (ii)  is  prototypical  of  a variety  of  inverse  wave  propagation  problems  of 
applied  mathematics.  We  refer  the  reader  to  (3),  [9],  for  instances. 

We  now  observe  that  problems  (i)  and  (ii)  are  equivalent.  In  fact,  if  we 
denote  by  L the  self-adjoint  ordinary  differential  operator  defined  by  the  boundary 
value  problem  (3.1)  (with  a boundary  condition  at  x = « supplied,  if  necessary),  then 
the  Riemann  function  is  just  the  distribution  kernel  of  cos  ti^L  , and  admits  the  spectral 
representation 

R(x,t,Xg  ,tg)  = / dp(X)  cos/\  (t-tQ)  $(x,\)  tUg,*). 

— CO 

In  particular 

£(t)  = R (0 ,t ,0 ,0 ) = / dp(X)  cos/x  t 
so  that  f is  the  Fourier  transform  of  p . Thus  knowledge  of  f and  knowledge  of  p 
are  equivalent,  so  problems  (i)  and  (ii)  are  equivalent. 

Since  problems  (i)  and  (ii)  are  equivalent,  the  solution  by  Gel'fand  and  Levitan 
(4)  of  problem  (i)  also  solves  problem  (ii).  We  prefer,  however,  to  deduce  the  solution 
of  problem  (i)  from  that  of  problem  (ii),  and  in  so  doing  present  a natural  interpreta- 
tion of  the  machinery  in  (4)  in  the  context  of  hyperbolic  p.d.e.  Besides  a better 
understanding  of  the  well-known  results  and  methods  of  (4] , we  are  immediately  led  to 
the  correct  stability  result  for  problem  (ii)  (Theorem  III) , and  to  the  solution  of  a 
number  of  other  inverse  boundary  value  problems  (see  [9)). 

We  should  point  out  that  the  equivalence  of  problems  (i)  and  (ii),  and  the  hyper- 
bolic interpretation  of  the  ideas  in  (4),  are  more-or-less  well  known.  However,  no 
careful  statements  on  the  lines  of  Theorem  II  have  appeared  in  the  literature,  nor  has 
a stability  result  of  the  type  of  Theorem  III  been  previously  asserted,  to  our  knowledge. 

Our  approach  to  problem  (ii ) is  based  on  several  elementary  properties  of  the  Riemann 
function,  especially: 


-7- 


(i)  (Group  property)  Let  U(t)  be  the  operator  which  maps  Cauchy  data 


f 


(u (x, tQ) (x»t^ ) ) for  a solution  of  (3.3)  at  time  tQ  to  the  Cauchy  data  for 
the  same  solution  (u(x,t+tQ) ,ufc (x,t+tQ ) ) at  time  t+tQ.  Then 

U(s)U(t)  - U (s+t ) (3.6) 

(ii)  (Progressing  Wave  Expansion)  The  distribution  R(y,t)  = R(y,t,0,0)  can 
be  decomposed: 

R(y,t)  = 6(y+t)  + 6 (y— t ) + K(y,t) 

I I 

where  K(y,t)  = K(y,-t)  has  one  more  derivative  than  the  coefficient  q in  the 
region  y <_  |t|,  vanishes  identically  outside  that  region,  and  on  the  boundary 
satisfies  the  transport  equation 

K(t,t)  - - j J*  q + h . (3.7) 

0 

These  assertions  will  be  made  precise  in  §4.  They  are  essentially  classical  results 
for  which  methods  of  proof  are  to  be  found  for  instance  in  Chs.  V and  VI  of  (1). 

Since  the  operator  U is  implemented  by  the  Riemann  function,  (i)  and  (ii) 
plus  some  symmetry  properties  of  R , together  imply  an  integral  equation  for  R : 

f (s , t ) = j-(f  (s+t)+f  (s-t))  = J°°  dy  R (y , t ) R (y , s ) . (3.8) 

2 0 

According  to  property  (ii ) , we  may  replace  the  upper  limit  of  integration  on  the  right- 
hand  side  of  (3.8)  by  max  ( |t  | , |s  | ) , and  we  obtain,  for  instance,  for  t >_  s > 0 

f(s,t)  = K (s , t ) + j*  dy  K(y,t)K(y,s)  . (3.9) 

0 

we  base  our  solution  of  problem  (ii)  on  this  nonlinear  Volterra  equation,  derived  in  a 
different  way  by  Gel'fand  and  Levitan.  The  hypotheses  of  Theorem  II  are  precisely  what 
is  necessary  to  ensure  that  (3.3)  has  a global  solution. 

Further,  our  method  of  solving  (3.9)  for  K is  manifestly  stable  against  small 
changes  in  the  data  f (Theorem  III).  Finally,  explicit  error  bounds  appear  (55)  which 
allcw  one  to  estimate  the  efficiency  of  our  approach  for  numerical  computation.  This 
circumstance  should  be  compared  with  most  other  treatments  of  problem  (i)  along  the  lines 
of  [4]  (see  especially  [6],  (31),  which  turn  on  the  solution  of  a linear  integral  equa- 
tion of  Fredholm  type  related  to  (3.9).  By  contrast,  our  error  estimates  involve  only  a 
l<*rer  bound  for  this  Fredholm  operator. 


Having  solved  (3.9)  for  K , it  is  relatively  easy  to  show  that  K is  the  appropriate 
piece  of  the  Rieraann  Function  for  a boundary  value  problem  (3.3),  with  coefficient  func- 
tion q related  to  K by  the  transport  equation  (3.7).  Ibis  is  accomplished  in  §7  . 
Finally,  Theorem  II  is  used  in  §8  to  supply  the  key  ingredients  of  the  proof  given  in 
16]  of  Theorem  I,  and  we  leave  the  matter  there. 


§4 . Properties  of  the  Riemann  Function 

This  section  is  devoted  to  the  necessity  part  of  Theorem  II,  that  is,  the  proof  that 
the  boundary  value  of  the  Riemann  function  of  a problem  of  form  (4.1)  below  must  indeed 
satisfy  the  conditions  (i)  and  (ii)  in  Theorem  II. 

Most  of  the  following  assertions  are  standard  and  can  be  found,  for  instance,  in  [1] , 
Chs.  V and  VI,  in  one  form  or  another.  In  a few  cases,  our  finite  diffe  entiability 
hypotheses  and  imposition  of  boundary  conditions  in  the  progressing  wav.’  construction  of 
R are  incompatible  with  readily  available  results;  however,  the  proofs  are  the  appropriate 
modifications  of  the  available  ones,  and  we  omit  them. 

Vfe  denote  by  ) the  solution  of  the  boundary  value  problem 


u.  - u + q u 
tt  xx 


u^  (0,t)  - h u (0,t) 


(4.1) 


with  initial  conditions 


R(X'WV  “ 4(X'V 


(4.2) 


3t  R(x-t'x0't0> 


0. 


t«t„ 


We  assume  that  q e w (R  ),  i > 1.  Then: 


I.  (Regularity)  Re  cf(R  x R+  x * ;8’  (R+  ))  n cm(R+  x x R i D' (1 R ) ) . 

t X0  t0  X X *0  V t 

Property  I follows  from  a theorem  on  solutions  of  (4.1)  with  smooth  Cauchy  data  and  the 
Schwarz  kernel  theorem,  via  standard  arguments. 
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II.  Symmetry) 


As  in  §3,  we  write  R for  the  distribution  kernel  on  R+  * R : 
R(y,t)  R(y,t;0,0) 


— (R(0,t-s)  + R:(0,t+s))  = / dy  R(y,s)  R(y,t) 

0 

Proof . By  III  , 


R(0,t+s)  = f dy  R(0,t+s,y,s)  R(y,s,0,C 


+ f“  dy  ( /t+S  do  R (0 , t+s  ,y , o ) ) D R(y,s,0,0).  (4.6) 

Os 

However,  using  the  symmetry  properties  (II, >, 

ft+S  do  R(0,t+s,y,o)  = J*  do  R(0,t,y,o) 
s 0 

D2R(y,s,0,0)  = -D  R(y,-s,0,0). 

Ulus  the  second  integral  (composition  of  kernels)  is  odd  in  the  parameter  s whereas 
the  first  is  even.  Replacing  s by  -s  in  (4.6),  adding,  and  multiplying  by  gives 


the  result. 


q.e . d. 


IV.  (Progressing  Wave  Expansion) 

R(y,t)  = 6 (y+t ) + 6 (y-t ) + K(y,t)  (4.7) 

where  K t (»'™+2'1(Q)  , 
loc 

Q = { (y , t ) e ]R+  x R ! y <_  |t|  ) , and  K(y,t)  = 0 , y > |t|  . 

Also  the  transport  equation  holds:  for  t i 0 , 

K(t,±t)  = - j J*  q + h . (4.8) 

0 

In  the  readily  accessible  literature,  the  proqressing  wave  expansion  is  usually 

CD 

developed  for  equations  with  c coefficients  and  no  boundary  conditions.  In  our  case, 
we  just  barely  get  away  with  it,  because  of  the  finite  differentiability  assumptions. 
Slightly  altering  the  notation  used  in  §3,  write 
R(0,t)  I 26  (t ) + ^(t) 

j (&(0,s+t)  + R(0,s-t))  » 6(s+t)  + 6(s-t)  + f(s,t) 
f(s,t)  - | $( s+t)  + hs-t))  . 

Then  ? c ' 1 ( R ) , f t ’ 1 (R  » R)  according  to  IV  ; according  to  II,  \ 

is  even  and  f is  symmetric. 
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In  view  of  all  this  symmetry,  we  consider  the  integral  relation  (4.5)  when  t, 

s ^ 0.  Then  the  progressing  wave  expansion  reads 

R(y,t)  = 6 (y-t)  + K (y,t) 

and  (4.5)  becomes , for  t > 0 , 

i(s-t)  + f(s,t)  = / dy(6(y-t)  + K (y,t))  (6  (y-s)  + K (y,s)) 

0 

= 6(s-t)  + K (s,t)  + f dy  K (y,t)  K (y,s) 

0 

or 

f(s,t)  * K (s,t)  + f dy  K (y , t ) K (y,s)  . (4.9) 

0 

Since  both  sides  of  (4.9)  are  continuous,  it  holds  for  t = 0 also.  Restricted  to 

0 < t < T , some  T > 0 , the  kernel  R defines  an  operator  of  the  form  I + KT,  where 

2 

1 is  the  identity  operator  on  L ([0,T]),  and  1^  is  the  integral  operator 

(y ) = /T  dt  K (y,t)  g(t) , 0 1 V L T • 

y 2 

Thus  I + K^,  is  a bounded  (Volterra)  operator  on  L ( [0 , T] ) , and  in  particular  is 
invertible. 

2 

Write  Ft  for  the  integral  operator  on  L ( [0 , T] ) with  kernel  f(s,t).  Wien  (4.9) 
can  be  written 

I + F “ (I  + K*)  (I  » 1^)  (4.10) 

which  shows  that  the  left-hand  side  is  a positive  semidefinite  symmetric  Fredholm  operator 
2 

on  L ( [0 , T]  ) . We  claim  that,  in  fact,  it  is  positive  definite.  Indeed,  it  is  the  product 
of  an  invertible  operator  and  its  adjoint,  and  therefore  positive. 

2 

Denote  by  e(T)  > 0 the  smallest  eigenvalue  of  I + F . Then  for  any  i|<  e L ([0,T]), 

( i(i,  (I+F  )i)/>  » /T/T  ds  dt  *(s)  i(i (t)  f(s,t)  + |i)i|2  ^ e(t)  f1  |i|i|2.  (4.11) 

L IM)  0 0 0 0 

This  completes  the  demonstration  of  the  necessity  of  hypotheses  (i)  and  (ii)  in 

Theorem  II.  The  next  three  sections  are  devoted  to  their  sufficiency. 
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§5.  Solution  of  a Nonlinear  Volterra  Equation. 

Let  T > 0 , and  set  = {(s,t)  :0^s^t^T}.  In  this  and  the  next  section, 
we  shall  show  that  the  equation 

f(s,t)  = K(s,t)  + Is  dy  K(y,s)  K(y,t)  (s,t)  e O (5.1) 

0 ~ 

has  a solution  K c (/*' 1 (o  ) for  every  symmetric  f e Jv'™'1  ( [0  ,T)  « [0 , T) ) which 

IOC  T IOC 

satisfies  condition  (4.11)  for  some  c (T)  > 0 . 

We  first  show  that  (5.1)  has  a continuous  solution  for  each  continuous  symmetric 
f satisfying  (4.11).  We  note  in  passing  that  continuous  solutions  are  necessarily  unique, 
since  (5.1)  is  a Volterra  equation. 

The  existence  proof  consists  of  three  steps. 

Step  1.  There  exists  some  t > 0 so  that  (5.1)  has  a continuous  solution  in 

This  is  true,  in  fact,  independently  of  the  hypothesis  (4.11).  We  denote  by  V the 

Volterra  operator  on  C°(Q^) 

V K(s,t)  = f(s,t)  - Is  dy  K(y,s)  K(y,t)  (s,t)  e Q-  . 

0 ~ 

Then  a standard  contraction  mapping  argument  shows  that,  provided  t is  small  enough, 
the  operator  V has  a fixed  point  in  C°(£h). 

Precisely,  we  obtain 

A.  Let  ||  ||-  denote  the  sup  norm  in  C°(Qj).  Then  the  ball 

Be,t(f)  = (g  e C°(Qj)  :||g-f||-  <.  e) 
is  invariant  under  V , provided 

2t  (e  + ||f  ||-)2  < e . 

B.  For  91i<32  « B^  j(f),  e as  above, 

II  % - v92  lit  ± 21 <E  + II f lit ) II  - g2n I • 

Thus  V is  a contraction  operator  for  t small  enough. 

C.  Denote  by  K the  fixed  point  of  V in  C°(Q^),  and  set 

K°  - f , K"  - ^ K° 

then 

||  K - I^||-  < {2t(e  + ||  f||- ))n  e (5.2) 

These  conclusions  are  completely  straightforward.  We  exhibit  them  only  to  make  explicit 
the  dependence  of  the  error  estimate  (5.2)  on  the  size  of  f . 


-13 


We  suppose  that  (5.1)  has  been  solved  in 


and  note  that  (5.1)  becomes,  in  the  region  0 < s < t <t<^T 
f(s,t)  = i(s,t)  + Js  dy  K (y , s ) 4(y,t) 


which  for  each  t c lt,T),  is  a linear  Volterra  equation  with  continuous  kernel,  hence 


has  a continuous  solution,  which  also  depends  continuously  on  the  parameter  t 


the  inhomogeneous  term  does 


Then,  in  the  notation  of  (4.10,11),  replacing  T by  t,  (5.3)  becomes 


and  (4.11),  implies 


We  observe  that  we  can  replace  e(t)  by  e(t'),  t'  ^ t,  in  (5.4),  in  particular  we  can 


replace  e(t)  by  e (T) . Together  with  the  Schwarz  inequality,  this  implies 


The  important  point,  of  course,  is  that  the  right-hand  side  of  (5.5)  is  independent  of  t 


Step  3.  For  t <_  s,t,  (5.1)  reads,  in  the  notation  of  Step  2 


This  problem  is  of  the  same  form  as  that  dealt  with  in  Step  1.  Therefore  we  can  solve 


(5.6)  in  a region  of  the  form  {(s,t)  :t^s^t^t+  5),  where  6 depends  as  in 


Step  1 on  the  uniform  norm  of  the  inhomogeneous  term 
f-(s,t)  « f(s,t)  - J*  dy  ♦ (y , s ) 4(y,t)  . 


Since  (5.5)  provides  us  with  an  estimate  of  the  size  of  f-  which  may  be  taken  inde 


indent  of  t , the  increment  6 may  also  be  chosen  independently  of  t 


Ifc  follo«s  that  a finite  number  of  repetitions  of  Steps  1-3  suffice  to  obtain  a 
. continuous  solution  K of  (5.1)  in  Q(T).  The  construction  proceeds  by  alternately 

solving  nonlinear  Volterra  equations  over  small  intervals,  and  linear  Volterra  equations 
. over  large  intervals.  The  error  in  each  step,  hence  the  cumulative  error,  and  the 

number  of  steps  required,  may  be  estimated  in  terms  of  the  three  numbers 

| l|f|lT,  A(T)  , e (T)  . 

Suppose  now  that  f e C "'(Q^).  By  differentiating  (5.1)  formally,  we  obtain  equations 
for  the  derivatives  of  K , which  we  solve  in  exactly  the  same  way  as  we  solved  (5.1). 

Thus  K e Cm(QT)  . 

We  remark  with  an  eye  to  the  proof  of  Theorem  III  that  the  following  easy  estimate 
shows  that  the  condition  (4.11)  is  stable  against  C°-small  perturbations  of  £ : 

I JT  ds  /T  dt  *(s)  (t ) (f  (s,t)  - f * (s  , t ) ) | 

0 0 

1 II  *11  % II  II  Hs,*)  - f (s,0  ||  II  2 

L [0,T]  I [0,T]  IT  [0 , T) 

s 

1 II  *H1 2,  T ||f  - f II  (5.7) 

L 10,T]  C (Qj,) 

for  any  ij/  e (0 , T) . 

It  follows  that,  for  any  f e which  satisfies  condition  (4.11),  there 

exists  a neighborhood  of  f in  Cm(QT)  in  which  condition  (4.11)  is  uniformly 
satisfied  (for  some  smaller  e > 0 , perhaps).  Easy  perturbation  arguments  show 
that  the  solution  K e C (Qp)  of  (5.1)  depends  continuously  on  the  left-hand  side  f, 
as  f ranges  over  any  such  neighborhood. 
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§6.  Existence  and  Stability  in 

The  existence  and  stability  of  solutions  of  (5.1)  with  locally  integrable  derivatives 
follow  from  the  basic  existence  and  stability  theorem  for  linear  Volterra  equations.  We 
state  this  elementary  result  for  the  reader's  convenience,  though  we  have  used  it  tacitly 
several  times  in  the  previous  section. 

Theorem . Suppose  a < b , Q = {(y,s)  : a ^ y < s j^b},  E a Banach  Space 

g e [a,b]  ; E) , V e Cm(Q;H).  Then  there  exists  a unique  if  e Cm([a,b);E)  with 

g(s)  = •fi(s)  + Is  dy  V(y,s)  v (y)  . 
a 

Furthermore,  the  solution  map: 

(f([a,b];E)  x C^(Q;K)  - cm(  [a,b]  ,E) 

(g,v)  m-  v> 

is  continuous. 

Finally,  there  exists  a W e Cm(Q;F),  called  the  resolvent  kernel  corresponding  to 

V , in  terms  of  which  the  solution  <f  may  be  written 

tf(s)  » g (s ) + Is  dy  W(y,s)  g(y)  . 
a 

Now  suppose  that  f eV^'1(B2)  satisfies  (4.11)  (The  case  f €h'"',1(B2)  , m > 2, 
loc  loc 

can  be  handled  ir.  exactly  the  same  way  so  we  give  details  only  for  m = 2.)  Since 

2,11  1 
W c c , the  previous  section  guarantees  that  the  solution  K of  (5.1)  is  in  C (QT> 

for  every  T > 0 , and  depends  continuously  on  f e C1 (QT> . 

2 

We  shall  first  shw  that  K has  a second  derivative  D^K  with  respect  to  its 

second  argument  in  C°([0,t):  i/dt.T]))  for  every  t e (0,T).  We  differentiate  (5.1) 

formally  twice  with  respect  to  t , and  obtain 

D2f(s,t)  - D2  K(s,t)  + J*  dy  K(y,s)  D2  K(y,t)  . (6.1) 

0 

New  f efc^^dl2)  means  that  D2f  € C°((0,tlgi  L*([t,T))t)  , for  each  t c (0,T). 

- 1 - 2 

Let  g:  [0,t]  -*•  L It,T)  by  g(s)  = D^ffs,*).  Then  the  basic  existence  theorem  asserts  that 
g(s)  « u(s)  t f dy  K(y,t)  u(y)  (6.2) 

0 o - 

has  a unique  solution  u e C ((0,t);  L^tt.T)).  It  remains  to  show  that  u may  be  iden- 
2 

tified  with  K as  a distribution.  In  terms  of  the  resolvent  Kernel  G corresponding 
to  K the  solution  u to  (6.2)  may  be  expressed: 
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u(s)  = g(s)  + J2 * * * * * 8  dy  G(y,s)  g(y)  . 

0 

The  kernel  G is  also  continuous.  If  we  regard  u and  g as  £>•  ( [t,T] )- valued 

functions  of  s e (0,t)  , then  for  0 « CQ((t,T))  we  may  write 

<u(s),0>  « <g(s),$>  + J8  dy  G(y,s)  <g(y),0> 

0 

- <f(s,*),0">  + j8  dy  G(y,s)  <f(s,*)0"> 

0 

■ < K(s,  •)  ,0") 

where  we  have  used  the  equation,  equivalent  to  (5.1): 

K(s,t)  « f(s,t)  + J8  dy  G(y,s)  f(y,t). 

0 

Thus  u is  indeed  the  second  derivative  of  K , and  we  conclude  that 
D^K  e C°([0,t]  i L1  [t , T] ) 

for  each  t e tO,T),  T > 0.  Finally,  since  K e C1  (QT>  depends  continuously  on 

2,1  2 
f € W (fij,),  the  stability  part  of  the  basic  theorem  implies  that  D^K  also  depends 

continuously  on  f , in  the  sense  of  C°([0,tl:  L^tt.T]),  for  each  t . 

The  other  two  formal  derivatives  of  (5.1)  read: 

D^f(s,t)  * D^K(s,t)  + ( x-K(s,s))  K(s,t) 

11  ds 

+ K(s,s)  Dx  K(s,t)  + D2  K (s  , s ) K (s , t ) (6. 

+ J8  dy  D^K(y,s)  K(y,t) 

0 

D D f(s,t)  = D D K(s,t)  + K (s  ,s  )D  K(s  , t ) 

11  (6. 

+ J8  dy  D K(y,s)  D K(y,t) 

0 ' 

From  these  equations,  which  hold  a priori  in  the  sense  of  distributions,  we  can 

2 i 

immediately  deduce  that  K e (/^^(Q))  we  leave  the  details  to  the  reader. 

2 

From  the  stability  assertion  about  D^K  and  some  obvious  fiddling  with  equations 

(6.3)  and  (6.4),  we  can  conclude  that  the  map  f K is  continuous  in  k'J^(Q),  in  a 
neighborhood  of  every  f satisfying  (4.11). 

Finally,  we  point  out  that  K inherits  from  f a property  somewhat  more  special 


than  that  of  being  in  fc'  , if  f is  of  the  special  form 

f(s,t)  - j <*(s+t)  + £(s-t)),  £ e W^(R>  • 

Indeed,  it  follows  from  a remark  in  §2  that  for  any  (a,b,c,d)  « R , the  function 
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§7.  Recovery  of  the  Boundary  Value  Problem 


We  shall  show  that  the  kernel  K constructed  in  the  last  section  is  the  regular 

part  of  the  Riemann  function  of  some  boundary  value  problem  (4.1),  thus  completing  the 

proofs  of  Theorems  II  and  III. 

Suppose  that  i1  £ w!"+^'*(R)  is  even.  Set 
loc 

f (s  , t ) - j(ha+t)  + £(s-t)> 

2(0)  - h . 

Then  f e is  syn,metric-  Suppose  that  f satisfies  (4.11)  for  each  T > 0. 

Let  K be  the  solution  of  (5.1)  constructed  in  §55,6,  with  f as  above.  Then  K e 


(Q)  , where  Q = { (s,t)  : 0 <_  s t). 
As  remarked  in  §6,  the  functions 


m+1 ,1 
are  in  w. 


s K K(as  + b,  cs  + d) 

of  their  appropriate  intervals  of  definition.  In  particular,  the  function 
q(s)  - -2  ± K(s,s) 


is  in  v/J1' 1 (»+  ) . 
loc 

We  claim  K solves  the  characteristic  mixed  problem 


{ “V  - -2-=-  + q(s)  > K(s,t 

L 3t  3s  ) 


s,t  £ Q 


(0,t)  h K(0,t)  » 0 (7.1) 

dS 

K(s,s)  = h - j Is  q . 

0 

Before  showing  that  K solves  (7.1),  we  remark  that  the  problem  (7.1)  has  a unique 
solution  in  ^loc'1  (Q)  , as  is  shown  by  standard  arguments.  Since  the  results  of  §4 
imply  that  the  regular  part  of  the  Riemann  function  of  the  boundary  value  problem  (4.1), 
with  q as  above , must  solve  (7.1)  also,  we  conclude  that  K is  the  regular  part  of  the 
Riemann  function.  We  have  therefore  proven  Theorem  II  if  we  shew  that  K is  a solution 
of  (7.1). 

To  do  this,  note  first  that  f solves 

7T  - 7T  f -°  ' f 5 0 ' 4 £ *• 


(Note  that  f has  second  derivatives  which  are  locally  inte9rable  in  * ) . According 
to  the  integral  equation 


we  have 
(1) 

(2) 


(3) 


f (s  ft)  - K(s,t)  + Is  dy  K (y , s ) K(y,t) 

0 


f(0,0)  =■  K(OfO)  = h 


1^-  (0,t)  - (0,t)  + K(0 ,0)  K(0  ,t) 

dS  dS 

* (0,t)  ♦ h K (0  , t ) :0  , t > 0 

dS 


3t2  3s2 


f (s,t) 


3t2  3s2 


K(s,t) 


S K(s,s)l 


K (s , t ) + 


- K (s , s ) D K(s  ,t ) - D K(s,s)  K (s  ,t)  - f dy  D2K(y,s)  K(y,t) 

r\  ^ 


(7.2) 


+ Is  dy  K (y , s ) D^K(y ,t) 
0 


The  differentiation  under  the  integral  sign  above,  and  various  integrations  by  parts  below 
are  justified  (for  the  only  questionable  case,  m = 1)  by 
Lemma 


Suppose  u c C°[t  ,t2l,  v = Du  e L^tt^ftj] 


Then 


/ 2 v « u(t2>  - u(t^) 


We  apply  this  Lemma  to  compute 
ay  r2 

0 


f3  dy  D K(y,s)  K(y,t)  » 

0 1 

I3  dy  D (D  K(y,s)  K(y.t)  - f dy  D^Jy.s)  D1K(y,t) 
ft  y ft 


* DjKfo'S)  K (s  , t ) - D K(0,s)  K (0  ,t ) 

+ fd  y K(y,s)  D2K(y,t)  - f dy  D (K  y,s)  D K(y,t)) 

0 1 0 y 
- D^K(s  ,s)  K (s  ,t ) - K (0  , s ) K(0,t)  - K(s,s)  D^s.t) 

+ K(0  ,s ) D K(0,t)  + /®  dy  K(y,s)  D2K(y,t)  . 

0 

Hence  for  0 < s <_  t 

K(s,s)  Dj^Kts.t)  - DjMs.s)  K (s  ,t ) 

» J8  dy  K (y ,s ) D2K(y,t)  - J8  dy  D2K(y,s)  K(y,t) 

0 1 0 

where  we  have  made  use  of  the  boundary  condition  (2).  By  virtue  of  the  identity 


(7.4) 
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equation  of  Volterra  type,  which  finishes  the  proof  of  (7.1)  in  that  case.  If  m = 0 , 
we  must  be  slightly  more  careful,  since  we  know  a priori  only  that  u is  locally  inte- 
grate in  Q = { (s  ,t ) : 0 <_  s t } . 

We  reason  as  follows.  Suppose  that  f is  approximated  by  f ( y 3' 1 . 

loc  n loc 

Then,  according  to  §§5,6,  the  corresponding  solutions  K e y3,1(o)  converqe  in 

n loc 

*l'  <Q>  to  K . As  noted  in  §2,  the  topology  of  y™'1  is  stronger  than  that  of  w"'1  . 

loc  loc 

If  we  denote  by  un  the  expression  (7.5)  formed  with  Kn  instead  of  K , the  sequence 
un  thus  converges  in  L\oc (Q)  to  u . However,  the  un  all  satisfy  (7.6)  and  are  con- 
tinuous, hence  vanish  identically.  Thus  u vanishes  almost  everywhere,  and  the  differ- 
ential equation  in  (7.1)  is  satisfied  in  the  sense  of  distributions. 

This  concludes  the  proof  of  sufficiency,  hence  the  entire  proof,  of  Theorem  II. 

To  prove  Theorem  III,  we  note  that  the  coefficient  q is  determined  by  the  transport 
equation  (characteristic  boundary  condition  in  (7.1)): 
q(s)  - -2  K(s,s). 

Since  it  was  shown  in  §6  that  K « V ' 1 (Q)  depends  continuously  on  h wm+1'1(R+) 

AOC  loc 

it  follows  that  the  map 

v/*+1,1((0,T]  ) 3 h h-  q < W^UO.Tl) 

is  (well-defined  and)  continuous  in  a neighborhood  of  every  satisfying  (4.11),  which 
is  exactly  the  assertion  of  Theorem  III. 
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§8.  Proof  of  Theorem  I 

Suppose  p : R -*  R is  a nondecreasing  function  satisfying  conditions  (i)  and  (ii) 
of  Theorem  1^ 

Lemma.  f , defined  as  in  Theorem  I,  cond—  (i),  satisfies  cond—  (ii)  of  Theorem  II. 

Proof.  Suppose  u e L2R  + has  support  in  (0,T).  Then 

/ dp(A)|u(A)|2  = lim  dp(A)  |u(A)  |2 

R N “ -N 

= lim  /N  d(  | A)  |u(X)  t2  + /N  do(A)  |u(A)  |2 
N-*od  -0  -N 


= f | u | 2 + lim  /N  do(A)(  JT  ds  u (s)cos  As)  (/Tdt  u(t)  cosAt) 

0 N-*“°  -NO  0 

where  we  have  used  the  ordinary  Parseval  formula  for  the  (cosine)  Fourier  transform. 

Now  according  to  condition  (i)  of  Theorem  I, 

lim  /N  dc(X)  cosA  s cosA't 
N -N 

exists  in  C°((0,T)  * (0,T) ) . Hence 

u(s)u(t)  f°  do(A)  cosAscosAt  = u(s)u(t)  f (s,t) 

- OD 

is  in  L'((0,TJ  * I0,T)),  and  by  Fubini's  Theorem 

T 

f“  do  ( A)  |u  ( A)  | = //  dsdt  u (s  )u  (t ) f (s  ,t ) . 

0 

Thus 

0 < f dr(A)|u(A)|2  » /T  I u 1 2 + /T/T  dsdt  u (s  )u  (t)  f (s  ,t ) 

” — 0 0 0 

and  condition  (ii)  says  that  the  inequality  is  strict  if  u is  not  identically  zero 
a.e.  This  however  is  equivalent  to  Condition  (ii)  of  Theorem  II,  as  remarked  in  §2. 

q.e.d 


We  may  therefore  apply  Theorem  II  to  conclude  the  existence  of  a hyperbolic 
boundary  value  problem 


u - U + qu  - 0 
tt  xx 


u^fO.t)  ♦ h u(0,t)  - 0 


(8.1) 


for  which  £ is  the  boundary  value  of  the  Riemann  function.  Here  q e (*+) 
Denote  by  $(x,A)  the  solution  (in  Cm+1(R+),  entire  in  A < C)  of 
-♦"(x,\)  + (q (x)  - \)$(x,A)  10  , x _>  0 
*(0,X)  « 1 , (0,X)  - -h  . 


(8.2) 
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Che  immediately  verifies  that  the  function 

u' (x,t)  5 cos/x  t 4,(x,  \ ) 

solves  (8.1)  with  initial  values 

u' (x, 0)  = $(x, X) 

X, 


It  follows  that 


In  particular 


ut(x,0)  ^ 0 


u (x,t)  = / dxQ  R(x,t ,x0,0) $(xQ, X)  . 


u A (0 , t ) = cos  A t 


with  K ( 

loc  * 


= / dxQ  R(0,t,xQ,0)  <MxQ,X) 

= / dx  R(x  ,t)  (|>(x  , X) 


= 4>(t,A)  + dxQ  K(xQ,t)  iHx0,X) 


As  in  §6,  denote  by  G the  resolvent  kernel  for  K . Then 

<Mt,X)  = cos/I  t + /*■  dx  G(x,t)  cos/x  x . (8.3) 

0 

the  other  hand,  K satisfies  equation  (5.1).  These  two  equations  ((5.1)  and  (8.3)) 
constitute  the  input  of  the  argument  in  16] , and  we  refer  the  reader  to  that  excellent 
reference  for  the  remainder  of  the  proof  of  Theorem  I. 

We  remark  that  equation  (8.3)  was  first  derived  by  Povzner  [7]  in  exactly  this  way. 


.Acknowledgement:  I am  pleased  to  thank  J.  Nohel  and  R.  Turner  for  several  helpful 

conversations . 
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